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Abstract 

We consider three-loop radiative-recoil corrections to hyperfine splitting in muonium generated 
by the diagrams with electron and muon vacuum polarizations. We calculate single-logarithmic and 
nonlogarithmic contributions of order a 3 (m/M)Ep generated by gauge invariant sets of diagrams 
with electron and muon polarization insertions in the electron and muon factors. Combining these 
corrections with the older results we obtain total contribution to hyperfine splitting generated by 
all diagrams with electron and muon polarization loops. Calculation of this contribution completes 
an important stage in the implementation of the program of reduction of the theoretical uncertainty 
of hyperfine splitting below 10 Hz. The new results improve the theory of hyperfine splitting, and 
affect the value of the electron-muon mass ratio extracted from experimental data on the muonium 
hyperfine splitting. 
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Muonium is a purely electrodynamic bound state, and existence of highly accurate ex- 
perimental results [1] makes it the best system for comparison of the precise quantum elec- 
trodynamic theory of hyperfine splitting with experiment 1 . The hyperfine splitting interval 
is proportional to the electron-muon mass ratio, and the current theoretical prediction is 

AE$ F8 (Mu) = 4 463 302 904 (518) (30) (70) Hz, (1) 

where the first error is due to the experimental error of direct measurement of electron-muon 
mass ratio m/M, the second error is due to the experimental uncertainty of the fine structure 
constant a, and the third error is an estimate of yet uncalculated theoretical corrections (for 
more details see [3] and [1]). The uncertainty of the electron-muon mass ratio dominates in 
the balance of errors, and therefore measurement of hyperfine splitting is the best source for 
the precise value of this mass ratio. We see from Eq. ([I]) that calculation of all theoretical 
corrections with magnitude above 10 Hz would improve accuracy of the electron-muon mass 
ratio. There are three series of yet uncalculated corrections of such magnitude pT]: a) 
single-logarithmic and nonlogarithmic radiative-recoil corrections of order a 3 (m/M)Ep, b) 
nonlogarithmic contributions of order (Za) 3 (m/M)EF, and c) nonlogarithmic contributions 
of order a(Za) 2 (m/ 'M)Ep [Z is the nucleus charge, Z = 1 for muonium). In this paper 
we complete the calculation of all single-logarithmic and nonlogarithmic radiative-recoil 
corrections connected with the electron and muon polarization loops. 

The radiative-recoil corrections of order a 3 (m/M)EF are enhanced by large logarithm of 
the muon-electron mass ratio M/m [5]. The leading logarithm cubed and logarithm squared 
contributions are generated by the graphs with the electron closed loops in Figs. [T]|4] (and by 
the diagrams with the crossed exchanged photon lines), and were calculated long time ago 
E]- Single-logarithmic and nonlogarithmic terms of order a 3 (m/M)Ep are generated by 
all diagrams in Figs. [TJjl], by the respective graphs with the muon loops, by the graphs with 
polarization and radiative photon insertions in the muon line, and also by the three-loop 
graphs with radiative photons in the electron and/or muon lines but without polarization 
loops. Below we calculate three-loop single-logarithmic and nonlogarithmic radiative-recoil 
corrections generated by the diagrams in Figs. [5j [6] both with the electron and muon loops. 
We combine these corrections with the earlier results, and present complete results for all 
radiative-recoil corrections generated by the diagrams with electron and muon polarizations. 




FIG. 1: Graphs with two one- loop polarization insertions 




See also recent discussion [2] of experimental feasibility of precise spectroscopy with another purely elec- 
trodynamic bound system, dimuonium. 
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FIG. 2: Graphs with two-loop polarization insertions 




FIG. 3: Graphs with radiative photon insertions 
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FIG. 4: Graphs with light by light scattering insertions 

There are numerous technical problems connected with calculation of the diagrams in 
Figs. [5j [6j As usual individual diagrams are ultraviolet and infrared divergent. In addition, 
the integrands in the Feynman parametrization for the exchanged integrals do not admit 
expansion in the small mass ratio before integration generating spurious divergences. This 
problem is solved by direct integration over exchanged momenta in four-dimensional Eu- 
clidean space. After angular integrations we used the method of overlapping small and large 
momentum expansions of the integrand to obtain systematic series over powers of the small 
electron-muon mass ratio and its logarithm [7J |SJ. Ultraviolet divergences are really not a 
problem and were dealt with by usual subtractions. 

By far the largest challenge is connected with the spurious infrared divergences. Gauge 
invariant sets of diagrams in Figs. [5j [6] are infrared finite, but each individual diagram 
is infrared divergent even in the infrared soft Yennie gauge. In order to obtain accurate 
results we need to separate the would be infrared divergent terms in the integrands, and 
explicitly cancel them analytically in the sum of contributions of different diagrams. In 
order to explain how such cancelation is achieved we start with the skeleton diagrams with 
two-photon exchanges in Fig. [7j All three-loop diagrams in Figs. |5j [6] may be interpreted 
as radiative corrections to the skeleton diagrams. We consider the momentum integrand 
in the skeleton product of the skeleton electron (muon) factor L^k) and the 

remaining part of the diagram. The factor L fJ/U (k) is the Compton scattering amplitude for 
virtual photons. Contributions to HFS generated by the diagrams in Figs. |5j|6] are described 
by the integrals similar to the integrals for the diagrams in Fig. [7j The only difference is 
that for the diagrams in Figs. [5j [6] we include in the integral radiatively corrected virtual 
Compton scattering amplitude instead of the skeleton one. A generalized low energy theorem 
holds for the virtual Compton scattering amplitude with subtracted anomalous magnetic 
moment contribution (see more on this subtraction below) [TJIH]. According to this theorem 
the electron (muon) factor L^ u {k) is suppressed by an additional factor k 2 /m 2 {k 2 /M 2 ) 
in comparison with the respective skeleton factor. Due to this suppression the integrals 
corresponding to the diagrams in Figs. [5| [6] are infrared finite, and the integration region 
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with momenta less than the electron mass are additionally suppressed. As a result, we can 
omit the atomic scale external virtualities of order ma, and calculate matrix elements in the 
scattering regime between the free electron and muon spinors. We use the Feynman gauge 
to obtain matrix elements of the gauge invariant sets of diagrams in Figs. [5j [6] We have 
derived some useful identities for the integrands that allowed to cancel the would be infrared 
divergent terms in the integrands before integration. Insertions of polarization operators in 
the diagrams in Figs. [5j [6] is taken care of by using the massive photon propagator for 
radiative photons (but not for exchanged photons) with the photon mass squared A 2 = 



4m 2 / (1 - v 2 ) or A 2 = 4M 2 /(1 



for the the electron and muon polarization loops, 



respectively. These massive propagators require an additional integration over velocity v 
with the weight J* dvv 2 (l — v 2 /3)/(l — v 2 ). 

The diagrams in Fig. [5] with the electron polarization loops generate nonrecoil and loga- 
rithm squared, single- logarithmic, and nonlogarithmic radiative- recoil contributions to HFS. 
It turns out that the gauge invariant anomalous magnetic moment in these diagrams does 
not generate radiative-recoil corrections (see, e.g., [HI E]). Then the radiatively corrected 
electron factor provides an additional suppression factor k 2 /m 2 in the skeleton inte- 
gral over the exchanged momenta. The wide integration region between the electron and 
muon masses m < k < M remains logarithmic. We calculated the nonrecoil contribution 
numerically, the logarithm squared and single-logarithmic terms analytically, and the non- 
logarithmic term numerically. The logarithm squared terms is already well known [BJ, and 
the single-logarithmic and nonlogarithmic contributions are as follows 



AE 



71 



53 
~6~ 



, M 

In — + 7.0807 

m 



a 3 m 



Ep. 



(2) 



The electron factor with muon polarization insertions in the diagrams in Fig. [5] provides 
an additional suppression factor k 2 /M 2 , and lifts characteristic integration momenta to the 
scale of the muon mass. Then these diagrams do not generate nonrecoil and logarithmic 
contributions to HFS. The respective leading recoil correction is a pure number, that we 
calculated it numerically 
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FIG. 5: Graphs with polarization insertions in the electron factor 

Consider now diagrams in Fig. [6j with radiative corrections in the muon line. The ra- 
diatively corrected muon factor provides an additional suppression factor k 2 /M 2 and the 
integral over exchanged momenta is nonlogarithmic. There is, however, another source 
of large logarithms in the case of electron polarization insertions in Fig. |6j The electron 
polarization insertions in the diagrams in Fig. [6j enter in the asymptotic regime, and the 
respective leading contribution to HFS is the product of the leading asymptotic term in 
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FIG. 6: Graphs with polarization insertions in the muon factor 
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FIG. 7: Diagrams with two- photon exchanges 

the high momentum expansion of the electron polarization operator and the radiative-recoil 
correction to HFS generated by the one-loop muon factor without polarization insertions. 
In other words the leading logarithmic contribution to HFS generated by the diagrams in 
Fig. [6] is obtained from the respective nonlogarithmic contribution of the diagrams without 
polarization insertions by substitution of the running coupling constant a(M) for radiative 
photons. We calculated also the nonlogarithmic contribution and obtained 
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The muon factor with muon polarization loops in the diagrams in Fig. [6] again provides 
the suppression factor k 2 /M 2 , but no enhancements. As a result these diagrams generate 
only nonlogarithmic radiative-recoil corrections. After numerical calculations we obtained 
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Other single- logarithmic and nonlogarithmic three-loop radiative-recoil corrections gener- 
ated by the diagrams with electron and muon polarization insertions were obtained earlier, 
and we collect respective results below. Single-logarithmic and nonlogarithmic radiative- 
recoil corrections generated by the diagrams with two electron or two muon loops in Fig. [I] 
are HDI 
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The diagrams with one electron and one muon loop in Fig. [8] produce only single- 
logarithmic and nonlogarithmic contributions to HFS [10J 
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The single-logarithmic and nonlogarithmic radiative-recoil contributions to HFS gener- 
ated by the diagrams in Fig. [2] with two-loop electron or two-loop muon polarization inser- 
tions were calculated analytically [10J 
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FIG. 8: Graphs with both the electron and muon (bold) loops 
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The single-logarithmic and nonlogarithmic radiative-recoil contributions generated by the 
diagrams with electron or muon polarization insertions in the exchanged photons in Fig. [3] 
have the form [11] 



AE 



a 3 m 



^ In - + 8.6945) ~E F . 
6 m J 7T" 3 M 



(9) 



+ 2 



+ 2 



FIG. 9: Muon line and electron vacuum polarization 



The diagrams in Fig. [9] with the electron or muon polarization insertions in the exchanged 
photons generate only single-logarithmic and nonlogarithmic radiative-recoil contributions 
to HFS [TTl 
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Combining all three-loop single- logarithmic and nonlogarithmic radiative-recoil correc- 



tions to hyperfme splitting due to electron and muon polarization loops in Eq. - Eq. (10) 
we obtain 
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For completeness let us mention that the three-loop radiative-recoil correction generated 
by the diagrams with one-loop fermion factors (and without polarization loops) in Fig. 10 
is also known [T2] 
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FIG. 10: Diagrams with two fermion factors 



The only still unknown single-logarithmic and nonlogarithmic three-loop radiative-recoil 
corrections are generated by the gauge invariant sets of diagrams with two-loop fermion 
factors without polarization insertions, and the diagrams with light-by-light insertions in 
the exchanged photons. Calculation of these corrections is a task for the future. 

The total result for all known three-loop single-logarithmic and nonlogarithmic radiative- 



recoil corrections to hyperfine splitting is given by the sum of the contributions in Eq. (11) 



and Eq. (12) 
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Numerically this contribution to HFS in muonium is 
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As was explained above the current goal in the theory of hyperfme splitting is to reduce 
the theoretical uncertainty below 10 Hz (see a more detailed discussion in [31 HU [13]). The 
muon and electron polarization operator contributions and other corrections collected in 
Eq. (13), together with the results of comparable magnitude in [HI [151 HB1 EH (HI EE] 
constitute a next step toward achievement of this goal. Phenomenologically, the improved 
accuracy of the theory of hyperfine splitting would lead to a reduction of the uncertainty 
of the value of the electron-muon mass ratio derived from the experimental data pQ on 
hyperfine splitting (see, e.g., reviews in [3"1 141 IT3]). 
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